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ABSTRACT 

Recent numerical simulations suggest that Population III (Pop III) stars were 
born with masses not larger than ~ 1OOM but typically ~ 4OM . By self- 
consistently considering the jet generation and propagation in the envelope of 
these low mass Pop III stars, we find that a Pop III blue super giant star has 
the possibility to raise a gamma-ray burst (GRB) even though it keeps a massive 
hydrogen envelope. We evaluate observational characters of Pop III GRBs and 
predict that Pop III GRBs have the duration of ~ 10 5 sec in the observer frame 
and the peak luminosity of ~ 5 x 10 50 erg sec -1 . Assuming that the E p — L p (or 
E p — E lt i so ) correlation holds for Pop III GRBs, we find that the spectrum peak 
energy falls ~ a few keV (or ~ 100 keV) in the observer frame. We discuss the 
detectability of Pop III GRBs by future satellite missions such as EXIST and 
Lobster. If the E p — E ly i SO correlation holds, we have the possibility to detect Pop 
III GRBs at z ~ 9 as long duration X-ray rich GRBs by EXIST. On the other 
hand, if the E p — L p correlation holds, we have the possibility to detect Pop III 
GRBs up to z ~ 19 as long duration X-ray flashes by Lobster. 

Subject headings: gamma rays: bursts — gamma rays: observations — gamma 
rays: theory 
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Introduction 



Gamma ray bursts (GRB) are the brightest phenomena in the universe. Long-soft type 
GRBs are considered to originate from deaths of massive stars such as Wolf-Rayet (WR) stars 
f lHjorth Sz Bloom! 201 if) . The most widely accepted s cenario for long GRBs is the collapsar 



scenario f jWoosley 



1993 



MacFadyen fc Woosleylll999l ). In this model, after the gravitational 



collapse of a massive stellar core, a black hole and an accretion disk system is formed and it 
launches a relativistic jet by magnetic field or neutrino pair annihilation process. If the jet 
can break out the stellar envelope successfully, a GRB is raised by converting the jet kinetic 
energy into the radiation energy. 

Owing to their brightness and detections at high redshift universe, GRBs are expected 
to be one of the powerful tools to probe the early universe. The development of observational 
instruments and early follow-up systems enable us t o discover some high re dshift GRBs. The 
most distant one ever is GRB 09042 9B at z = 9.4 (Cucchiara et al.l 



at z = 8.3 (e.g. iTanvir et al.ll2009l ; ISalvaterra et al 



2009 



2011) and GRB 090423 



Chandra et al.l 120101 ) follows it. If 



GRBs can be raised by first stars and detected, we will draw informations about the early 
universe, e.g., the star formation history and the reionization history. 

First stars in the universe, so called Population III (Pop III) stars, are considered to 
be formed from metal free gas in the very early universe. The metal-free primordial gas 
cools less efficiently compared to the metal-contained present-day gas, which allows the pri- 
mordial gas to have larger fragmentation masses. Since it was considered that the whole 



(Abel et al. 


2002; 


Bromm et al. 


2002) 



dieted to be very massive ~ 100 — 1000i\ 
recent studies suggest that this is not always the case an d that a massive gas clump can ex- 



perience further f ra gmentation to f orm a bina ry system (jTurk et a. 



Clark et al.ll201lh . iTan fc McKed (|200J) and McKee & Tan! (12008 ) suggested that the UV 



20091 : IStacv et al.ll2010 



radiation from the central protostar can ionize the surrounding neutral gas and suppress 
the accretion onto the protostar. They analytically investigated this feed back effect on the 
protostar evol ution and found that P op III stars finally obtain mass typically ~ 14OM . 
More recently, iHosokawa et al.l (1201 ll ) performed two-dimensional simulations of the proto- 
star evolution including the above feed back effect and found that Pop III stars finally obtain 
masses typically ~ 4OM . They also concluded that the UV radiation from the central star 
eventually stops the mass accretion and the growth of the star by the evaporation of the 
surrounding gas. 

It is considered that metal free Pop III stars do not lose mass, keeping large hydrogen 



envel opes until the pre-supernova stage, because of the low opacity envelopes (fWooslev et al. 



20021 ). The final fate of a Pop III star depends on the stellar mass (IHeger et al.ll2003l ). After 
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the stellar core collapse, those Pop III stars in the range of 1OM < M < 25M would 
explode as supernovae and form neutron stars as remnants. Those in 25M < M < 4OM 
form black holes as remnants after the fall back accretion of the envelopes onto the temporally 
formed neutron stars. More massive stars (40M Q < M < 14OM and 260M o < M) would 
fail to blow out their envelopes and promtoly form massive black holes, except those stars 
in 140M Q < M < 26OM who end as pair-instability supernovae due to the explosive 



Fryer et al. 


2001; 


Suwa et al. 


2007) 



There have b een some studies about the productivity of GRBs from massive Pop HI 



stars (> lOOM-^ (iMeszaros fc Reed l20ld iKomissarov fc Barkovl 12010 : ISuwa fc Iokal 12011 
Nagakura et alJl2012h . The former two studies assumed a massive black hole surrounded by 



an accretion disk as an outcome of a massive stellar collapse, and estimated the accretion 
rate onto the black hole. Then they evaluated the jet luminosity and showed that the burst 
activity of a Pop III star is observable by current detectors. In ISuwa &: Iokal (1201 ll ). they 
analytically studied the jet propagation in the stellar envelope and showed that massive Pop 
III stars (~ 9OOM ) can produce GRBs although they have large hydrogen envelopes, since 
the long last ing accretion provides enough energy and time for the successful jet breakout. 
In addition, iNagakura et al.l (120121 ) performed two-dimensional relativistic hydrodynamic 
simulations in which the accretion onto a black hole and the jet production are treated 
in a self-consistent way for stellar models of massive Pop III stars (915M ), Wolf-Rayet 
stars (initially 16M Q ), and low mas s Pop HI stars (4OM ). They confirmed the validity of 
the analytic results in ISuwa &: Iokal (120 111 ) and also found that 4OM Pop III stars can be 
progenitors of GRBs, but did not study their observational characteristics and detectability. 



The idea of GRBs from blue super giants (BSG) was suggested in IMeszaros &: Rees 



(120011 ). Although they considered the jet dynamics in the stellar envelope, they treated a 
steady jet and did not reflect the central engine activity caused by the change of the accretion 
rate. They did not evaluate th e possibility of GRBs from BSGs quantitatively. On the other 
hand, IWoosley Sz Hegerl (120121 ) discussed gamma-ray transients from Pop III BSG collapsars 
by investigating the mass accretion of the outermost layers of a star, but did not discuss 
the jet propagation and the jet break out. Assuming that the conversion efficiency of the 
accretion energy to the radiation energy ~ 10~ 2 , they found that Pop III BSGs can produce 
long gamma-ray transients with duration 10 4-5 sec and luminosity io 48-49 erg sec -1 . In this 
paper, we simultaneously investigate both aspects (the jet propagation and the central engine 
activity) in a self-consistent way by including the following physical processes; the stellar 
collapse, the non-steady jet injection, and the jet propagation in the stellar envelope. By 
doing this, we quantitatively discuss the possibility of the jet break out and GRB especially 
for low mass Pop III stars (around 40 M ). 



-4- 



In §2, after introducing the stellar models and the jet propagation models, we investigate 
the productivity of a GRB focusing on a 4OM Pop III star, whi ch is a Pop HI star with the 
typical mass reported by the state-of-the-art simulation done in iHosokawa et al.l (120 111 ). In 
§3, we calculate the observational characters, such as the duration T go , the peak luminosity 
L p and the spectrum peak energy in the observer frame E° hs , of GRBs from 40M© Pop III 
progenitors. Then we evaluate the detectability of such Pop III GRBs by future detectors 
such as Lobster and EXIST in detail, varying the redshift of a burst. We apply the above 
discussions to different progenitor models with masses of 30 — 90M Q . In the last part of §3, 
we evaluate the light curves of Pop III GRB radio afterglow emissions and their detectability 
by the Low Frequency Array (LOFAR) and the Expanded Very Large Array (EVLA). §4 is 
devoted to the summary and discussions. 



2. GRBs from low mass Pop III stars 
2.1. Progenitor and Relativistic Jet models 



We employ a pre-collapse stellar model of z40.0 by lWoosley et al.l (120021 ). which provides 



the structure of a 4OM with zero metallicity at the final phase of the stellar evolution. It 



is considered 


that 


(Heeer et al. 


2003) 



Then we consider the subsequent evolution of the stellar collapse and 



the je t propagation in the stellar envelope following the similar prescription to ISuwa fc Ioka 



( 1201 it ) 



We first assume that the collapse proceeds in a spherically symmetric manner without 
pressure support so that each mass shell of the star within mass [M r , M r + dM r ] and radius 
[r, r + dr] falls into the central core in the free-fall time scale t^{r) = y/ r 3 /GM r . We 
calculate the mass accretion rate from the expression M(r) = dM r /dtg(r). When the mass 
of the central core becomes 3M , we identi fy that a black hole is formed since the maximum 
possible mass of a neutron star is ~ 3M Q (IRhoades fc Rumnilll974l ; IChitre fc Hartld 119761 ). 
After the formation of the black hole, we assume that a cold relativistic jet with a constant 
opening angle 9j = 5° has been launched and we take this moment as the origin of time 
(f = 0). 

There are mainly two candidates for the jet production mechan isms, i.e., the neutrino 
annihilation process and the magnetic process. In lSuwa fc Iokal (120111 ). they showed that the 
jet model based on the neutrino annihilation process is not appropriate for producing GRBs 
from very massive Pop III stars. Accordingly we adopt the jet model based on the magnetic 
process. In this model, the jet injection luminosity Lj Ct (t) is considered to be represented as 
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Ljet(t) = TjM{t)c 2 (see ISuwa fc Iokal l201ll and references there in), where the constant rj is 



an energy conversion efficiency and we take the value of t] = 6.2 x 10 -4 . This is a calibrated 
value so as for Wolf-Rayet stars to reproduce the energetics of canonical local long GRBs, 
i.e., nearly 10 52 ergs of energy should be injected into the relativistic jet after the breakout. 

In this paper it looks like that we neglect the effect of th e stellar rotat i on an d treat 
the stellar collapse in spherically symmetric way. According to iKumar et al.l (120081 ) , when 
the stellar rotation is taken into consideration, the accretion time scale onto the central 
BH for each mass shell t acc ( r ) can be represented as t acc ( r ) ~ tfi( r )/&, whe re a ~ 0.1 is 



the st andard dimensionless viscosity parameter of the disk. As described in ISuwa fc Ioka 
( 1201 ll ). we regard that this uncertain factor is absorbed within the calibrated parameter rj. 
Therefore, we think that the disk formation is implicitly taken into account and that the 
spherically symmetric prescription makes sense. 

In the following sections, we consider the propagation of a jet in the stationary stellar 
envelope. For the 4OM Pop III model, the He core mass Mn e and the He core radius rn c 
are Mh c ~ 22M Q and rn c ~ 10 11 cm, respectively. Therefore, the collapse time scale of the 
He core is estimated as £ co n ~ 600 sec. On the other hand, the time scale for the jet head to 
reach the outer edge of the He core can be evaluated as t cr0 ss ~ ^He/O-lc ~ 30 sec, since we 
can see that the average velocity of the jet head within the He core is ~ 0.1c (see Fig. [p. 
Accordingly, t cross <C t co ii holds. We confirm that this inequality holds better in outer layers 
and that calculations based on the stationary envelope is self consistent. 



2.2. Jet Propagation in the Pop III Star Envelope 

First, we consider the propagation of a jet in the stellar envelope. As pointed out in 
the previous subsection, we approximate that the stellar envelope is stationary and that the 
density profile is the same as that in the pre-supernova stage until the jet break out. A 
jet propagating through the stellar envelope forms forward and reverse shocks at its head. 
Here, we assume that the separation between these shocks is small compared to the distance 
from the stellar center. From the continuity of the momentum flux at the jet head, we have 



flMatznerl 120031 ) 



PjC %{Yf h y% - i3 h y + p 3 = P ^K(v h f3 h y + p„ (i) 

where p,h,T, (3 and P represent the density, the specific enthalpy, Lorentz factor, the velocity 
divided by the speed of light c, and the pressure, respectively. The subscripts j, h and * 
stand for the jet, the jet head and the stellar envelope, respectively. We consider a cold jet 
so that we can neglect Pj in l.h.s. of Eq. ([!]). In r.h.s. of Eq. ([!]), we can approximate 
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h* ~ 1 and neglect P*, since stellar material is non-relativistic. Then the velocity of the jet 
head ((3 h ) is expressed as 



7rr fe 2 g/p,(r fe )c 3 
[Ljit-rh/iPjc)) 



1/21 - 1 



(2) 



where Lj(t — r^j {(3jc)) and refer to the jet luminosity and the position of the jet head, re- 
spectively. We use the formula Lj = Ti{rh9j) 2 pjC 2 hjY 2 f3jC in calculating Eq. ([2]). Accordingly, 
the position of the jet head is calculated as r h (t) = £P h (t')cdt. 

The jet head consists of shocked stellar matter and shocked jet material. They are 
relativistically hot and expand sideways of the jet forming a cocoon. We assume that almost 
all the jet energy goes through the shocked region into the cocoon during the jet propagating 
in the stellar envelope. The cocoon expands laterally by balancing its pressure with the ram 
pressure of the stellar matter as 

P c = pAPl + P* ~ p*c 2 /3 2 , (3) 

where the subscript c refers to the cocoon and P c ^> P* is assumed. Since the cocoon 
consists of relativistically hot materials, P c can be expressed as P c = E C /(3V C ), using the 
cocoon volume V c and the cocoon energy E c . Now, we suppose the shape of the cocoon as 
a cone, then V c (t) = 7ir 2 (t)rh(t)/3. In our jet model, the cocoon energy can be expressed 
as E c (t) = r]M acc (t)c 2 , where M acc (t) is the mass accreted to the black hole by the time 
t. Substituting all these expressions into Eq. (J3]), the cocoon expansion velocity (3 C can be 
calculated as 

„,.-, r h (t) / 4»)M, cc (i) 

where M(r/ l ) = (47rrjj/3)p*. Then the position of the cocoon edge is given as r c (t) = 
$1 p c (t')cdt' '. 

Now we discuss whether Pop III stars can raise GRBs by following the time evolution of 
the positions of the jet head and the cocoon edge. If the jet head reaches the stellar surface 
earlier than the cocoon edge, we consider that the star can raise a GRB since we can expect 
a successful jet breakout. On the other hand, if the cocoon edge reaches the stellar surface 
earlier, we can expect that the mass accretion is suppressed and that the relativistic jet is 
stalled on the way. This looks like a failed GRB. 

Fig. [1] shows the time evolution of the jet head velocity (the red solid line) and the 
cocoon velocity j3 c (the green dashed line). In this figure, the time variability of velocities 
comes from the discontinuity of the stellar density profile and the mass accretion rate. As 
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pointed out in e.g.. IMeszaros fc Reesl ( 1200 ll ). we can see that the jet head accelerates drasti- 
cally after entering the hydrogen envelope. In addition, we find that the jet head propagates 
faster than the cocoon edge all the way through the stellar envelope except for the very early 
time. We also find that the jet head breaks out of the stellar envelope ~ 400 sec after the 
central engine is activated. Thus, we conclude that a 4OM Pop III star has the possibility 
to raise a GRB. 

Note here that 4OM Pop III stars are thought to end their lives as blue super giants 
(BSG), keeping large hydrogen envelopes with radii ~ 10 12 cm. This i s because the opacit y 
is too low to induce the mass loss from metal free stellar envelopes flWoosley et all 120021 ) . 
In general, the progenitor of a local long GRB is not considered to be a super giant star 
with a hydrogen or helium envelope but to be a Wolf-Rayet star with no hydrogen or helium 
envelope and radius ~ 10 10 cm. The observational reason is that every supernovae associating 
with long GRBs belongs to type Ibc. Theoretically, it is considered tha t a super giant star 
has a too largely extended envelope for the jet to break out successfully (IMatznerl 120031 1 and 



it cannot raise a GRB. From the results here, however, we confirm that a BSG 
is compact enough for a successful jet breakout. 



10 



12 



cm 



We should note that IMeszaros fc Reesl (120011 ) and IWoosley fc Hegerl (120121 ) have sug- 
gested the possibility of GRBs or gamma-ray transients from BSGs. However the former 
treated a steady jet and did not quantitatively evaluate the possibility, while the latter fo- 
cused on the formation of the accretion disk around the BH and did not discuss the jet 
break out. In this section we quantitatively confirm the possibility of GRBs from BSGs by 
consistently considering the stellar collapse, non-steady jet injection, and the jet propagation. 



3. Observational Properties of Pop III GRBs 

3.1. The prompt emission 

In this section, we consider observational characters of Pop III GRBs. We assume that 
soon after the jet breakout, the jet emission can be seen as a GRB and the burst lasts 
until the whole stellar envelope accretes completely. We suppose that the efficiency for 
converting the jet energy to the radiation energy is 10 %. Accordingly, we can calculate 
expected properties of the burst, such as the peak luminosity (L p ), duration (T 90 ) and the 
isotropic energy (i? 7iiso ). Note here that we estimate the T 90 as the period during which 90 
% of the burst's energy is emitted. Furthermore, we evaluate the time-integrated spectrum 
using empirical laws for GRBs. There are some correlations which hold between the time- 
integrated spectral peak energy in the observer frame (E° hs ) and the peak luminosity (L p ) 
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or the isotropic energy (-B 7i i S o)- One is t he E p — L p correla tion ( lYonetoku et al.ll2004f ) and 

iso correla 

two correlations are represented as 



the other is the E p — E ltiso correlation fjAmati et al.l 120021 ) . The functional forms of these 



10 52 erg sec 



-l 



2 x 10 -5 



£° bs ( 



1 + z) 



1 keV 



2.0 



E° p bs {l + z) 
1 keV 



80 



E 



7, iso 



10 52 erg 



0.57 



(5) 



(6) 



respectively. 



Table [T] shows the results of our model for a 40M Q Pop III star. As can be seen from 
Table [TJ Pop III GRBs radiate as much energy (E 7 ^ ao ~ 10 54 erg) as the most energetic local 
long GRBs do, while L„ ~ 5 x 10 50 erg sec -1 is smaller by a factor of ~ 10. Moreover, the 



duration (~ 10 5 sec) of Pop III GRBs is much longer than that of local long GRBsjj All 
these differences come from the fact that although the progenitor of a local long GRB has no 
hydrogen envelope and is more compact, a Pop III progenitor has a large hydrogen envelope. 
Since a Pop III star experiences no mass loss and keeps a more massive hydrogen envelope, 
the energy supply to the central engine can last longer time. This enables the central engine 
to be kept active for much longer time. This also enables the burst to have a much longer 
duration and to have the vast isotropic energy. On the other hand, it takes longer time for 
the jet to break out the larger stellar envelope and the jet energy after the breakout is more 
damped. This causes the Pop III GRBs to have lower luminosities. 

We evaluate the observed peak energy for either the case that the E v — L p correlation 
holds or that the E p — E ltiso correlation does. For the case of E p — L p correlation, E° hs ~ 5 
keV is in the X-ray region, whereas for the E p — E^^ so correlation, the peak energy in the GRB 
frame is larger than that of a local long GRB, because of the larger -E^iso value. However, 
the cosmological redshift effect reduces the peak down to E° hs ~ 120 keV, which is slightly 
softer than that of a local long GRB. 



1 In lSuwa fc Iokal (|201l[ ). they evaluated the duration of the Pop III GRB from a 915M Q star as ~ 1, 500 
sec in the GRB frame, but this is a wrong value. We find that the correct value is ~ 15, 000 sec in the GRB 
frame, which is similar to the one obtained here (~ 6000 sec in the GRB frame). This is because the larger 
mass of the 915A/© Pop III star compensates with its larger radius (~ 10 13 cm), as we can see from the 
expression of the free-fall time. 
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3.2. The detectability of Pop III GRBs 



In this subsection, we discuss the detectability of Pop III GRBs. In lSuwa fc Iokal ( 1201 ll ). 
they found that Pop III GRBs are too dim to trigger Swift Burst Alert Telescope (BAT). We 
obtain the similar conclusion even we employ the different progenitor from them. Therefore, 
we here discuss whe ther Pop III GRBs trigger the future satellite missions such as Lobster 
rtGehrels et allkoii ) and EXIST^ in detail. While Lobster will have energy window range 



of 0.3 - 5 keV, EXIST will have that of 5 - 600 keV. 

An event is regarded to be detected if the number of the signal photons within the 
detector energy range [E min , E max ] satisfies the following relation, 

C +A '^KJ<^ > (S/N)min (7) 

Here, to is the time when an event comes in the detector's field of view and the detector 
starts to observe the event, At is the exposure time for the event, and t ODS is the time 
from the beginning of the burst in the observer frame. iV S i g , A and N\, g refer to the signal 
photon number flux, the area of the detector and the number flux of background photons, 
respectively. (S/N) min is the critical signal to noise ratio needed for detection. Assuming 
that the background photon flux is constant and using the signal photon energy flux within 
the detector energy window range [E^, E max ], Eq. ([7]) can be written as 

/ sig (t ,At) >/ sen (At), (8) 

where ^ ^ 

/ sig (t , At) = — (yj 

is the energy flux averaged over the exposure time and 

rt +At , , , 

/sen(At) - J ; o+At 7 s ° bs (S/N) min A-^< 2 At-^ (10) 

is the energy flux sensitivity within [E min , E max \. In our model, / S i g (^obs) can be calculated 
from 

( . L^)S B E ::EN { E )dE 2:cc _, 

4^4 j" °° EN(E)dE g ' [ ' 



2 http://exist. gsfc.nasa.gov/ 
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In Eq. ffTTj) . L 7iiso (i bs) is the isotropic equivalent luminosity of the burst at t bs- di is the 
luminosity distance calculated with cosmological parameters (f2 m , Q\) = ( 0.28, 0.72) and the 
Hubble parameter H = 70 km sec -1 Mpc" 1 . N(E) is the Band spectrum feand et al. 1993 ) 



with the typical parameter values, a = — 1 and = —2.3. We discuss the detectability of 
Pop III GRBs for either the case that the E p — L p correlation holds or that the E p — E^^ ao 
correlation does. 

First, we consider the case of the E p — L p correlation. Since E° bs ~ 5 keV, in this 



case L obster is more appropriate for detecting Pop III GRBs. Recently, iGhirlanda et al. 



( 120101 ) studied the time dependent spectral characteristics of several individual bright GRBs. 
They found that the isotropic equivalent luminosity L 7i i SO (/obs) correlates with the time 
resolved spectrum peak energy -Ep(t b s ) for each GRB and that the functional form of the 
correlation is very similar to the time integrated E p — L p correlation (Eq. 03)). Note that 
they calculated the time resolved spectrum by integrating the signal flux within 1 sec time 
bin around each time. Accordingly, if we assume the validity of the time- resolved E p (t h s ) — 
-^7,iso(^obs) correlation, which is obtained by replacing L p and E° hs in Eq. ([5]) by L 7 4 SO (£obs) 
and E° hs (t b s ), we can discuss the detectability using the condition in Eq. OH]). 

The Lobster sensitivity for a soft source (a power-law photon index of —2) is estimated 
to be 1.3 x 10~ n erg c m" 2 sec" 1 (0.3 — 5 k eV, 5a) at one calendar day (an effective exposure 



time of ~ 2500 sec; see lGehrels et al.ll2012l ). On the other hand, the sensitivity for a proposed 
exposure time per pointing in a realistic operation (~ 450 sec) is calculated to be 3.1 x 10~ n 
erg cm -2 sec -1 (0.3 — 5 keV, 5<r). The assumed spectral parameter in this estimation is 
reasonable for a GRB with E° hs ~ a few keV. We discuss the detectability of a Pop III GRB 
by Lobster using the sensitivity in At ~ 450 sec as a realistic case and in At ~ 2500 sec as 
an optimistic case. 

In Fig. [21 we compare the energy flux of GRBs from 4OM Pop III stars with the 
detection thresholds of Lobster. The abscissa is the time from the beginning of a GRB, i.e. 
from the jet break out, in the observer frame. The green, sky-blue and the blue solid lines 
represent / S i g (t b s ) of Pop III GRBs at z — 9, 14 and 19, respectively, calculated from Eq. 
QTT]) . The red and magenta dashed lines correspond to / sen of Lobster in a realistic case 
(magenta) and an optimistic case (red). From Fig. [2J we can see that / S i g (^obs) does not 
change significantly over At ~ 450 sec or ~ 2500 sec around each time, so we can approximate 
/sig ~ const, over the considered exposure times. Then, Eq. ([8]) can be rewritten as 

/rig(*0) > /scn(At). (12) 

Eq. ffl2|) indicates that if /sig(^o) when an event comes into the Lobster field of view is larger 
than / sen (At) for given At, we can observe the event from to to to + At. From Fig. [21 we 
find that Pop III GRBs at z — 9, 14 and even at z — 19 have the possibility to trigger 
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Lobster. Lobster will detect a Pop III GRB as a long duration X-ray flash with nearly 
constant luminosity 

Subsequently, we consider the case of the E p — E~ /tiso correlation. Since E° hs ~ 120 keV, 
in this case EXIST is the better instrument for detection. Note that the E p —E 1 ^ so correlation 
is the correlation between the total radiated energy and the time-integrated spectrum, we can 
regard E° hs as the observed peak energy time- averaged within each burst. So, we evaluate 
/sig (5 — 600 keV) assuming that the spectrum is the Band type with E° hs ~ 120 keV, a = — 1 
and (3 = —2.3 and that the spectrum does not change with time. The sensitivity of EXIST 
for a proposed exposure time in the longest time-scale at the on-board process (At ~ 512 



sec) is calculated to be f scn ~ 2.4 x 10 10 erg cm 2 sec 1 (5 — 600 keV, 5a) flHong et al. 



2009|). 



We show the results for 4OM Pop III stars in Fig. El Here again, the abscissa is the time 
from the beginning of a burst, i.e. from the jet break out, in the observer frame. The green, 
sky-blue and the blue solid lines represent / s i g (t bs) of Pop III GRBs at z = 9, 14 and 19, 
respectively, calculated from Eq. (TIT]) . The red dashed line represents the EXIST sensitivity 
described above. In this case also, / S i g (t bs) is approximately constant over At ~ 512 sec, so 
Eq. (jSJ) can be rewritten in the form Eq. f fl2|) . From Fig. [31 we can see that although Pop 
III GRBs at z = 14 and 19 do not trigger EXIST, Pop III GRBs at z = 9 have the possibility 
to trigger EXIST. EXIST will detect such a Pop III GRB as a long duration X-ray rich GRB 
with nearly constant luminosity. 



3.3. Other progenitor models 



McKee Sz Tan! (120081 ) and iHosokawa et all (120111 ) studied the mass of a Pop III star at 



its birth by calculating the evolution of a primordial protostar in analytical or numerical way. 
Because the initial angular momentum of a primordial gas cloud is considered to be large 
enough, in the star formation phase, a protostar and a circumstellar accretion disk system 
is formed and the protostar gains mass by the accretion of the surrounding gas through 
the disk. They found that the UV radiation from the protostar eventually stops the mass 
accretion and the growth of the star by evaporating the surrounding gas and eventually the 
disk. They also found that the final mass of a Pop III star depends on the degree of the 
angular momentum transport within th e accretion disk and on the magnitude of the initial 
angular momentum of a gas cloud. In IHosokawa et al.l ( 120111 ). the degree of the angular 
momentum transport is characterized by the ap-par ameter of the disk, whe re the larger ao 
value means the larger mass accretion. In fig. SI of IHosokawa et al.l (1201 ll ). a Pop III star 
finally gains ~ 50M Q with a = 1.0, ~ 40M Q with a — 0.6 (fiducial case), and ~ 35M Q 
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with cto = 0.3, for the fiducial magnitude of the initial angular momentum. Furthermore, 
the mass of a Pop III star depends on the initial angular momentum of the star-forming gas 
cloud and it gets ~ 85M when the initial angular momentum is reduced by 30 % of the 
fiducial one with a = 0.6 (fiducial case). 

Accordingly, in this subsection, we investigate whether 30 — 9OM Pop HI stars can be 
the progenitors of GRBs. The stel lar models of 30 — 40M^ a re given by lWoosley et al.l (120021 ) 
and those of 41 — 9OM a re fromlrleger fc Woosley ( 2010 ). In Fig. HJ we show the density 
profiles of selected models. IWoosley et al.l ( 120021 ) showed that all the 30 — 40M^ P op III stars 
end their lives as blue super giants (BSG). According to lHeger fc Woosleyl ( 120101 ) . although 
41 - 44, 60 and 70M o Pop III stars end as BSGs 45, 50, 55, 65, 75, 80, 85 a nd 90M© Pop III 
stars end as red super giants (RSG). As shown in lHeger fc Woosleyl (120101 ). the RSG branch 
in the higher mass stars appears due to the primary nitrogen production in the hydrogen 
burning shell. 

For these stellar models, we investigate whether these stars can raise GRBs by consid- 
ering the jet propagation in the stellar envelope with the entirely similar manner as in §2. 
Note here that although for a Pop III star with M < 40M Q the fall bac k effect should be 
taken into consideration in the formation of a black hole remnant (see e.g.. iMacFadyen et al. 



200ll ; iKumar et al.l 120081 ). we neglect this effect in this section and discuss it in §4. For 
BSGs (30 — 44, 60 and 70 M ), the jet head propagates faster than the cocoon edge almost 
all the way through the stellar envelope like in Fig. [TJ On the other hand, for RSGs, we find 
that the jet head reaches the surface as early as or even later than the cocoon edge does. 
Therefore, we conclude that although Pop III BSGs (30 — 44, 60, 7OM ) have the possibility 
to raise GRBs, Pop III RSGs do not. From the above discussions, we can say that although 
Pop III BSGs (~ 10 12 cm) are compact enough for successful jet breakouts, Pop III RSGs 
(~ 10 14 cm) have too largely extended low density envelopes for jets to break out successfully 

(Fig. m). 

Subsequently, we evaluate the observational characters and the detectability of these 
Pop III GRBs at z = 19. The results for some models are shown in Table [2j We can see the 
same characteristics for Pop III GRBs as described in §3.1. Here again, we consider either 
the case that the E p — L p correlation holds or that the E p — i? 7i i so correlation holds. First, we 
consider the case of the E p — L p correlation. In Fig. the red and the magenta dashed lines 
represent .f SCTI (At) (0 .3-5 keV, 5a) of Lobster in the optimistic and realistic case, respectively 
(IGehrels et al.ll2012f ). The green, blue, sky-blue, grey and black solid lines correspond to 
the energy flux f sig (t ) of Pop III GRBs at z = 19 for 30,40,44,60 and 7OM progenitors, 
respectively. The abscissa is the time from the beginning of each burst, i.e. from the jet 
break out. We find that while it is difficult for Pop III GRBs from 60 and 70 M progenitors 
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to trigger Lobster, Pop III GRBs from < 44M stars have the possibility to trigger Lobster. 
Second, we consider the case of the E p — E ltiso correlation and Fig. [6] shows the results. The 
red dashed line represents the EXIST sensitivity. Here also, the green, blue, sky-blue, grey 
and black solid lines correspond to the energy flux f S i g (t ) of Pop III GRBs for 30, 40, 44, 60 
and 70 M Q progenitors, respectively, but at z = 9. We find that only Pop III GRBs from 
< 44M Q stars have the possibility to trigger EXIST. 

From Fig. [5] and Fig. [6J we can see that the energy fluxes of Pop III GRBs from < 44M 
progenitors are larger by a factor of 3 to 4 than those from 60, 70M Q stars. This is because 
more massive progenitors exhibit larger radii (see Fig. H|). The larger radius a progenitor 
has, the longer time it takes for the jet to reach the stellar surface and the more damped the 
jet luminosity after the breakout is. 



3.4. The afterglow 



In this subsection, let us discuss the afterglow of a Pop III GRB following iToma et al. 



(120 111 ). In the external shock model of an afterglow, we consider that a relativistic ejecta 
with isotropic equivalent kinetic energy Ei SO and a half opening angle 6j moves through the 
interstellar medium with density n making a shocked region at the head of it. In the shocked 
region, some fractions, e# and e e , of the internal energy are provided to the magnetic field 
energy and the energy for the electron acceleration, respectively. The accelerated electrons 
are assumed to have a number distribution of the form N(-j e ) oc 7~ p and to raise afterglow 
emissions through the synchrotron radiation and the inverse Compton emission. As fiducial 
parameter values of the external shock model, we adopt -Ei so ~ 10 55 erg, n = 1 cm~ 3 , 
e e = 0.1, = 0.01, 9j = 0.1 and p = 2.3. In §3, we assumed that 10 % of the jet energy 
was converted into the energy for the prompt photon emission. Accordingly, we consider 
that the remaining 90 % of the jet energy is used for afterglow emissions. As we saw in §3, 
E y jsn ~ 10 54 erg, so here we adopt E iso ~ 10 55 erg as a fiducial value. We refer to lToma et al. 



(120 111 ) for other parameter values. Under this model, we calculate the afterglow light-curves 
at 10 GHz, 1 GHz and 100 MHz. We show the results in Fig. [7] in the case of a Pop III 
GRB at z — 19. The red, green and the blue solid lines refer to the light curves at 10 GHz, 
1 GHz and 100 MHz, respectively. 

Let us discuss the detectability of such Pop III GRB radio afterglow emissions. The Low 
Frequency Array (lofar|| has frequency coverage from 10 to 250 MHz and the detection 
threshold of 0.2 mJy (la level) at 100 MHz for 1 hr integration time. From Fig. 0, we find 



3 http://www. astron.nl/ 
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that the 100 MHz radio afterglow emission is not detectable by LOFAR. On the other hand, 
the Expanded Very Large Array (EVLA) has frequency coverage from 1 to 50 GHz and the 
detection thresholds of 5. 5 p,Jy and 1.8 /xJy (la level) at 1 GHz and 10 GHz, respectively 
for 1 hr integration time ( jPerley et al.l 1201 if ). From Fig. we can see that the energy fluxes 
at 1 GHz and 10 GHz are much larger than the detection thresholds of EVLA. Once a Pop 
III afterglow emerges, we have the possibility to detect it at any time by EVLA, even if it 
occurs at such a high redshift universe (z = 19). 



4. Summary and discussion 



GRBs are the brightest phenomena in the universe. Their detections at high z universe 
[z ~ 9) motivate us to expect GRBs to be one of the powerful tools to probe the early 
universe. Focusing on the high z universe, we s hould consider the ass ociation of GRBs with 
Pop III stars. Recent numerical simulations (IHosokawa et al.l 120111 ) suggest that Pop III 
stars obtain mass typically ~ 4OM at their birth. Zero metallicity stars are considered 



not t o lose mass during entire life because of the low opacity envelopes ( jWoosley et al. 



20021 ). Theref ore, they enter i nto th e pre- supernova stage k e eping large hydrogen envelopes. 
According to IWoosley et al.l ( 120021 ) and iHeger fc Woosleyl ( 120101 ) . Pop III stars end their 
lives as BSGs or RSGs depending on the amount of primary nitrogen produced in the shell 
burning. 

In this paper, we investigate whether such low mass Pop III stars ranging from 30 to 
9OM can be progenitors of GRBs. For this purpose, we consider the jet propagation in the 
stellar envelope and analytically calculate the evolution of the jet-cocoon structure. In BSG 
envelopes, the jet head velocity is larger than the cocoon velocity all the way except for the 
very early time and we can expect a successful jet breakout. On the other hand, in RSG 
envelopes, the cocoon edge reaches the stellar surface as early as or even earlier than the jet 
head. We confirm that Pop III RSGs have enough la rgely extended e nvelopes for jets to be 
stalled on the way and cannot raise GRBs as shown in iMatznerl (120031 ) . We also confirm that 
Pop III BSGs are compact enough for the successfu l jet breakout and have the po ssibility to 
raise GRBs as suggested in iMeszaros fc Reesl ( 1200 ll ) and lWoosley fc Hegerl (120121 ) . It should 
be noted that the BSG models from Woosley and H e ger u sed above ignored the effect of 
the rotation on the stellar evolution. lEkstrom et al.l ( 120081 ) found that when the rotation 
is included, Pop HI s tars within our noticed mass range end up as RSGs not BSGs. But 
Ekstrom et al.l (120081 ) considered the evolution of a star with an extrem ely high rotation 



velo city as one half th e critical velocity. Recent cosmological simulations ( IStacy et al.ll2010 



and I Clark et al .1120 111 ), however, suggested that Pop III stars are born in binary systems. In 
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this case, the angular momentum which the star forming gas clump initially has is divided 
into the spin of each star and the orbital angular momentum so that these stars may rotate 
less rapidly. Therefore, we think that such rapidly rotating stars they considered are rare 
and the calculations based on BSG models from Woosley and Heger make sense. 

Using our model, we evaluate observational characters of Pop III GRBs. We predict 
that although Pop III GRBs radiate as much energy as the most energetic local long GRBs, 
Pop III GRBs are slightly less luminous than local long GRBs due to their much longer 
burst duration. Assuming that the E p — L p (or E p — E 1)iso ) correlation holds for Pop III 
GRBs, we predict that Pop III GRBs have the much sof t er (or mildly softer) spectra than 



local long GRBs in the observer frame. IWoosley fc Hegerl (120121 ) predicted that the gamma- 



ray transients from low metallicity BSGs have duration of 10 4-5 sec and the luminosity of 
10 48-49 erg sec -1 , similar to Pop III GRBs considered here. Their transients are fed by 
the mass accretion of the outer most layers of stars and the accretion rate onto the BH 
(~ 1O~ 4 M sec -1 ) is much smaller than that in the Pop III GRB case, which imply typically 
from - 1O" 3 M sec 1 to ~ 10 2 M & sec l . Moreover, they simply assumed a central engine 
model with the roughly estimated conversion efficiency from the mass accretion to the jet 
energy as ~ 0.1, while we consider a central engine which is driven by the magnetic process 
implicitly taking the disk accretion into account with the efficiency of rj ~ 6.2 x 10" 4 @ 
Note that we choose this value so as for Wolf-Rayet stars to reproduce the energetics of 
local long GRBs. Therefore, although t he characters are s imilar among them, we expect 



that gamma-ray transients considered in IWoosley fc Hegerl ( 120121 ) are different events from 
GRBs considered in this paper. 

We also discuss the detectability of Pop III GRBs by future satellite missions such as 
Lobster and EXIST in detail. If the E p — E ltiso correlation holds, we have the possibility to 
detect Pop III GRBs at redshifts z ~ 9 as long duration X-ray rich GRBs by EXIST. On 
the other hand, if the E p — L p correlation holds, we have the possibility to detect Pop III 
GRBs up to z ~ 19 as long duration X-ray flashes by Lobster. 

We b riefly comment the ex pected observable GRB rate per year by Lobster using the 
results of de Souza et al. (12011 ) . We calculate the observed GRB rate per year dN G ^ B /dz 
as 

d^GKB _ ^obs dN GRB 
dz Air dz 

where dN GRB / dz, fl ohs and ?7beam correspond to the intrinsic GRB rate (the number of on-axis 
and off- axis GRBs) per year, the detector field of view and the beaming factor of the burst. 



4 Note that IWoosley fc Hegerl (|2012T ) considered the rotationally supported disk structure and the mass 
accretion from it so that the meaning of the conversion efficiency is different from ours. 
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In Fig. 6 of Ide Souza et al.l (1201 if ), they showed gL/Vgrb/g^ for an optimistic case and we use 

ilso adopt the values of r/beam ~ 0.01 and f2 D bs ~ 0.5 sr for Lobster 



their values. 1 


iere, a 


( Gehrels et al. 


2012) 



and 0.4 Pop III GRBs per year at z — 9, 14 and 19, respectively. 

At last, we briefly discuss employed assumptions in this paper. Firstly, we assume that 
all the stars considered in this paper (30 — 9OM Pop III stars) form black holes directly 
after the stellar core collapse (see §3.3). It should be noted, however, that this is not always 
the case especially for less massive stars. Shortly after the onset of the core collapse, a 
neutron star and a shock wave, which propagates outward or is stalled, are considered to 
be formed at first. Behind the shock wave, a fall back accretion of the shocked envelope 
could be present and the continuous accretion onto the neutron star eventually leads to a 
black hole formation. Although the early activity of the central engine could be affected by 
the accretion details, i.e., the direct accretion or fall-back accretion, the conclusion of this 
paper is hardly changed. This is because the mass accretion at the interested time in this 
paper is coming from the massive envelope so that the central engine already collapsed to a 
black hole at the corresponding time. In addition, since the energy budget of the shock head 
is dominated by the envelope accretion, the details of the early phase does not affect the 
shock evolution in the late phase. For a more massive (> 40M Q ) star, on the other hand, 
the energy of the shock wave is too low to explode even the portion of the envelope, so a 
black hole would be formed directly and our assumption is fully justified in this case. Note 
that the mass threshold bet ween the di rect or fall-back induced black hole formation is still 
under the debate (see e.g. lFryerlll999l ) and beyond the scope of this paper. Secondly, we 
assume that the whole stellar envelope accretes onto the BH (see §3.1). In order to confirm 
the validity of this assumption, we evaluate the binding energy of each layer of the stellar 
envelope and compare it with the typical energy of a supernova outgoing shock ~ 10 51 erg, 
which is injected around ~ 10 km from the center. We find that the binding energy becomes 
larger than 10 51 erg within r < 5 x 10 9 cm. This means that the outgoing shock s hould stall 
on the way and we expect little mass ejection. Recently, iQuataert fc Shiodd (120121 ) suggested 
that in the late stage of the stellar evolution, the pre-supernova burning leads to a significant 
mass ejection from the outer envelope. But they considered only the case of a 4OM star 
with metallicity Z = 10 -4 and commented that the amount of mass ejection depends on the 
metallicity and rotation etc. Thus, the amount of the ejectable mass is uncertain for the 
progenitors employed here so that we do not consider this effect in this paper. 
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Table 1: The comparison of typical long GRBs and Pop III GRBs with 40 M ( 





long GRB 


Pop III GRB (z = 9) 


-E 7 ,iso [erg] 


1Q52-54 


10 54 


L p [erg sec^ 1 ] 


^q52-53 


6 x 10 50 


T 90 [sec] 


10 1 " 3 


6 x 10 4 


E° bs [keV] 


1Q 2-3 


5.5 (E p - L p ) 
120 {E p — E 1)iso ) 



Table 2: The observational characteristics of Pop III GRBs at z = 19 for various progenitor 
masses. 



mass [Mq] 


30 44 60 70 


E 7iiso (10 54 ) [erg] 
L p (10 50 ) [erg sec" 1 ] 

T 90 (10 5 ) [sec] 
E° h * [keV] (E p - L p ) 
E; bs [keV] (E p - £ 7 , iso ) 


0.94 1.1 1.5 1.6 
5.1 6.4 1.7 1.9 
0.87 1.3 10 11 
2.5 2.8 1.4 1.5 
54 60 70 75 
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time [sec] 

Fig. 1. — The time evolution of the jet head velocity flh (the red solid line) and the cocoon 
velocity j3 c (the green dashed line) during the propagation in the stellar envelope. The 
abscissa is the time from which the jet is activated. The progenitor is a 4OM Pop III star. 
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Fig. 2. — The comparison of the energy flux / s i g (to) f° r a P°P 111 GRB with the detection 
thresholds of Lobster, f SCQ (At) (0.3-5 keV, 5a). The abscissa is the time from the beginning of 
a GRB, i.e., from the jet break out, in the observer frame, to is the time when the event comes 
into the Lobster field of view and Lobster starts to observe the event. At is the proposed 
exposure time of Lobster. The green, sky-blue and the blue solid lines correspond to / sig (0.3- 
5 keV) for the GRB from a 40M Q Pop III star at z = 9, 14 and 19, resp ectively. The red an d 
the magenta dashed lines represent / sen (At) (0.3-5 keV, 5cr) of Lobster (IGehrels et al .1120121 ). 
corresponding to the optimistic case of At ~ 2500 sec and the realistic case of At ~ 450 sec, 
respectively. If / S i g (t ) > f scn (At) holds, we regard that Lobster observes the Pop III GRB 
from to to to + At. 



-23 - 



10" 



10 



■10 



> 

CD 

o 
o 

CO 

uS 1 







■11 



10 



•12 




10 £ 



10' 



10 4 
time [sec] 



10' 



10 c 



Fig. 3. — Same as Fig. [2] but for EXIST case. EXIST will have the limited en- 
ergy range of 5-600 keV. The red dashed line represents the EXIST sensitivity / sen ~ 
2.4 x 10 -10 erg cm -2 sec -1 (5 -600 keV, 5a) in the longest exposure time-scale at the on- 



board process (At ~ 512 sec) (IHong et al.ll2009l ). 
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Fig. 4. — The density profiles for some of the Pop III stellar models (jWoosley et al.l I2002 



Heger fc Woosleyll2010l ). Those who end their life as BSGs correspond to 3OM (red), 40M Q 
(green), 44M Q (blue), 6OM (magenta) and 7OM (sky-blue), respectively. They have radii 
±Qi2-i3 cm. Those who end their life as RSGs correspond to 45M (grey) and 65M (black), 
respectively. They have radii ~ 10 14 cm. We can see that RSGs have much more extended 
and lower density envelopes than BSGs have. 
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Fig. 5. — The comparison of the energy flux f S i g {to) for Pop III GRBs with the detection 
thresholds of Lobster, f sen (At) (0.3-5 keV, 5a). The abscissa is the time from the beginning 
of a GRB, i.e., from the jet break out, in the observer frame, to is the time when the 
event comes into the Lobster field of view and Lobster starts to observe the event. At is the 
proposed exposure time. The green, blue, sky-blue, grey and black solid lines correspond 
to / sig (0.3-5 keV) of Pop III GRBs at z = 19 from 30,40,44,60 and 7OM progenitors, 
respectiv ely. The red and t he magenta dashed lines represent / sen (At) (0.3-5 keV, 5cr) of 
Lobster (IGehrels et al.ll2012l ). corresponding to the optimistic case of At ~ 2500 sec and 
the realistic case of At ~ 450 sec, respectively. If f S i g (t ) > / sen (At) holds, we regard that 
Lobster observes the Pop III GRB from t to t + At. 
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Fig. 6. — Same as Fig. [51 but for EXIST (5-600 keV) case. The red dashed line represents 
the EXIST sensitivity / scn ~ 2.4 x 10~ 10 erg cm -2 sec -1 (5 -600 keV, 5a) in the longest 
exposure time-scale at the on-board process (At ~ 512 sec) flHong et al.l 120091 ) . Note that 
we focus on Pop III GRBs at z = 9 in this figure. 
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Fig. 7.— The afterglow light-curves of a Pop III GRB at 100 MHz, 1 GHz and 10 GHz in 
the case that the GRB happens at z = 19. The red, green and the blue solid lines correspond 
to the light-curves at 10 GHz, 1 GHz and 100 MHz, respectively. The abscissa is the time 
from the beginning of the prompt emission, i.e. from the jet break out time, in the observer 
frame. Here, we adopt the parameter values of n — E 55 = e e -i = €3,-2 = f(p) = — 1 
and p = 2.3. Q x denotes Q/10 x . 



